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RELIABILITY  DATA  ANALYSIS  MODEL 


INTRODUCTION 


Current  Army  regulations  1 1  I  reouire  the  specification  of  Quanti¬ 
fied  reliability  goals  in  the  development  of  new  weapon  systems, 
fhese  reliability  goals  must  also  be  verified  bv  tests  prior  to  final 
acceptance  and  fielding  of  the  system.  Reliability,  bv  definition, 
is  a  probabilistic  quantity.  In  addition,  the  amount  of  testing  that 
can  be  conducted  is  limited  by  cost  and  time  considerations  and  conse¬ 
quently  reliability  assessment  must  be  conducted  within  a  statistical 
framework.  To  assist  in  this  assessment  effort  a  statistical  data 
analysis  computer  model  has  been  prepared  and  is  described  in  this 
report.  Essentially,  the  model  computes  Point  estimates  and  confi¬ 
dence  limits  for  reliability  of  components,  subsystems  and  system 
from  component  test  data. 

The  specific  reliability  index  considered  is  the  mission  relia¬ 
bility  defined  as  the  probability  that  a  given  system  will  success¬ 
fully  perform  its  intended  function  without  failure  for  a  specified 
period  of  time  under  a  given  set  nf  conditions.  Time  can  be  given 
in  terms  of  clock  time,  rounds,  cycles,  miles,  etc.  The  system  may  be 
required  to  perform  a  number  of  missions  over  its  expected  life. 

^Arrny  Regulation  AR70R-50,  "\rmv  Materiel  Reliability  and  Maintain- 
ability,"  Headquaters  Dept .  of  the  Army,  Washington,  D.C.,  8  January 
1968. 
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PURPOSE  OF  THE  RELIABILITY  COMPUTER  MODEL 


The  reliability  model  is  intended  to  serve  a  number  of  useful 
purposes.  First,  reliability  analyses  for  the  realistic  situations 
encountered  in  weapon  system  testing  and  for  all  but  the  simplest 
assumptions  made  for  computation  entail  rather  complicated  mathemati¬ 
cal  formulation  and  tedious  computation.  The  computer  model  conse¬ 
quently  provides  a  systematic  means  of  incorporating  more  complex 
analyses  with  less  restrictive  assumptions  into  a  routine  data  analy¬ 
sis  procedure.  Second,  the  model  provides  greater  speed  in  obtaining 
results  of  data  analyses.  Third,  the  mathematical  procedures  and  for¬ 
mulations  incorporated  into  the  model  provide  a  means  of  standardizing 
reliability  assessment  among  design  agencies,  test  agencies  and  the 
user,  finally,  the  model  provides  a  basis  for  determining  the  type 
and  amount  of  test  data  required  to  perform  a  particular  kind  of 
analysis.  This  provides  insight  into  planning  of  tests  and  data 
collection  procedures  required. 


TYPE  OF  DATA  CONSIDERED  FOR  ANALYSIS 


The  type  of  data  considered  for  analysis  is  failure  and  suspen¬ 
sion  data  on  a  component  or  a  number  of  components  making  up  a  sub¬ 
system  or  system.  Definition  of  failure  depends  on  the  particular 
requirements  and  the  mission  profile  of  a  given  svstem  and  can 
include  part  breakage,  out-of-tolerance  performance,  incipient 
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failure,  safety  hazard,  etc.  Suspended  data  is  obtained  when  a  compo¬ 
nent  is  removed  from  the  test  without  failure,  the  test  is  terminated 
without  failure  or  when  the  failure  of  a  comnon^nt  is  attributable  to 
the  failure  of  another  component  or  cause. 

Since  failure  rates  and  resulting  mission  reliabilities  are 
generally  transient,  all  data  is  assumed  collected  as  a  function  of 
system  age.  Sneci f i cal ly ,  the  age  given  in  terms  such  as  rounds,  time 
or  miles  on  individual  components  within  a  system  at  failure  or  sus¬ 
pension  is  required.  In  general  it  is  not  enough  to  record  onlv  the 
total  test  time  and  number  nr  fai lures  excent  in  the  case  of  constant 
failure  rate. 

TYPE  OP  INFORMATION  CENERATEn 

The  tvpe  of  information  generated  hv  the  model  includes  results 
of  analyses  for  each  individual  connon^nt  within  a  subsvstem,  for 
each  subsystem  making  up  a  given  system  and  ror  the  entire  system. 
Component  information  includes  results  of  distribution  selection 
procedures,  Kolmogorov-Smimov  goodness-of-fi t  tests,  maximum  likeli¬ 
hood  estimates  of  distribution  parameters,  noint  estimates  o^  the 
mission  reliability  as  a  function  of  svstem  age  ror  the  non-constant 
failure  rate  case  and  confidence  limits  on  average  mission  reliabilitv 
over  the  system  life  for  the  constant  failure  rate  case.  Subsystem 
and  system  information  generated  includes  r>rmot  estimates  and  confi¬ 
dence  intervals  on  mission  reliability  for  the  constant  failure  rate 
case. 
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A  more  complete  description  of  required  input  data  for  the  model 
and  the  resulting  output  information  will  be  presented  in  subseouent 
sections  dealing  with  the  computer  program. 


USUAL  ASSUMPTIONS  FOR  DATA  ANALYSIS  AND  THEIR  LIMITATIONS 

It  is  worthwhile  at  this  point  to  review  the  usual  assumptions 
made  in  performing  a  reliability  assessment  for  a  given  system  and 
the  desirability  of  relaxing  these  assumptions  for  some  of  the  realis¬ 
tic  cases  encountered  in  weapon  system  testing.  Undoubtedly,  the  most 
common  assumption  made  in  reliability  data  analysis  is  the  assumption 
of  constant  failure  rate  for  components  and  systems.  The  most  signi¬ 
ficant  implication  of  this  assumption  is  that  the  probability  of 
failure  of  a  component  or  system  is  independent  of  its  age.  Components 
are  assumed  to  fail  at  completelv  random  points  in  time.  This,  of 
course,  does  not  permit  the  full  treatment  of  the  cases  of  early 
failures  and  wear-out  failures  common  to  mechanical  components. 

There  are  a  number  of  reasons  why  the  constant  failure  rate 
assumption  is  so  often  made.  First,  the  constant  failure  rate  is 
often  assumed  primarily  to  simplify  analysis  and  computation.  Second, 
straightforward  techniques  for  determining  confidence  intervals  on 
component  reliability  from  test  data  are  readilv  available  for  con¬ 
stant  failure  rate.  Third,  testing  procedures  are  greatly  simplified 
since  only  the  total  number  of  failures  and  total  test  time  are  re¬ 
quired  for  data  analysis.  Also,  the  amount  of  data  required  for  each 
component  is  not  as  critical  as  in  those  cases  where  theoretical 
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distribution  selection  must  be  considered.  Final.lv,  for  systen  rolin- 
bilitv  verification  tests,  individual  component  data  is  not  required 
with  all  failures  being  treated  as  system  failures.  This  is  true 
onlv  if  all  components  have  equal  test  times. 

Assuming  a  constant  failure  rate  can  lead  to  three  types  of 
errors;  the  first  is  in  the  computation  of  reliability  for  fixed  'ITBr 
(mean  time  between  failures),  the  second  is  in  the  computation  of 
confidence  limits,  and  the  third  is  in  the  estimation  of  'TBn  when 
suspension  or  censored  data  has  been  generated.  The  magnitude  of  the 
error  in  computing  reliability  for  fixed  '1TBF  cannot  be  determined  ;n 
general.  For  illustrative  pur^ses,  however,  consider  a  hvpotbeticoi 
case  in  which  a  system  is  made  up  of  eoual  components  in  series. 

Assume  further  that  each  component  has  a  WcibuJ 1  distribution  and  an 
'1TBF  (mean  time  between  failure)  eoual  to  twice  the  expected  system 
life*  The  number  of  missions  over  system  life  is  assumed  to  be  150. 

The  magnitude  of  the  error  in  assuming  constant  failure  rate  can  now 
be  determined  for  this  particular  case*  Table  I  lists  the  results  of 
average  system  reliability  for  different  numbers  of  components  in  the 
system  and  for  various  values  of  the  Weibull  shape  parameter  R.  For 
a  shape  parameter  B  equal  to  1.0  the  Weibull  distribution  reduces  to 
the  exponential  which  describes  the  distribution  for  constant  failure 
rate.  For  values  of  R  greater  than  1.0  the  failure  rate  is  an  increas¬ 
ing  function  of  time  characteristic  of  wear-out  phenomena.  The  greater 
the  value  of  R  the  more  peaked  the  distribution  is  about  the  mean. 

As  can  be  seen  from  Table  I,  large  differences  can  be  obtained  when 
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AVERAGE  SYSTEM  RELIABILITY  FOR  EQUAL  COMPONENTS  HAVING  WEI  BULL  DISTRIBUTION  OP 
FAILURE  TIMES. 
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0.999993  0.999934  0.999672  0. 090345  0.996727 
0.999996  0.999962  0.999812  0.999625  0.998124 


assuming  exponential  or  constant  failure  rate  for  Weibull  components 
that  in  actuality  have  3  values  greater  than  1.0.  ^or  ex ample,  a  sys¬ 
tem  of  SO  equal  components,  each  component  having  a  Weibull  shape  para¬ 
meter  equal  to  3.0  would  have  an  average  system  reliability  of  0.072. 
Assuming  a  constant  failure  rate  op  the  other  hand  would  yield  a  com¬ 
puted  reliability  of  0.84b  which  is  significant ly  different  than  the 
true  value  for  this  case. 

Reliability  for  the  Weibull  distribution  withOM  was  computed 
using  renewal  theory  where,  in  this  rase,  a  component  is  replaced  or 
renewed  upon  failure.  The  mathematical  formulation  for  the  renewal 
case  is  described  in  the  Analytical  "ethods  section  of  this  report. 

The  second  source  of  error  in  assuming  constant  failure  rate  is  in 

the  computation  of  confidence  limits.  Consider  as  an  example  problem 
a  test  sample  oF  failure  times  given  as  2000^  2200,  1800,  1000  and 

2100.  For  this  problem  the  MTBF  is  2000  and  the  standard  deviation  is 
158,  Lower  confidence  limits  on  the  true  *TTBF  can  readily  be  computed 
assuming  different  underlying  distributions  of  failure  times.  ror 
example,  the  lower  00°/  confidenced  means  for  the  exponential  (constant 
failure  rate)and  the  normal  (increasing  failure  rate)  distributions 
are  1250  and  1758  respectively.  It  is  common  in  reliability  verifica¬ 
tion  tests  to  renuire  that  the  lownr  confidenced  '1TBF  (or  reliability) 
exceed  a  given  fixed  value  for  acceptance.  If  the  requirement  in  the 
above  problem  were  1500,  the  assumption  of  the  exponential  distribution 
would  lead  to  rejection  and  the  assumption  of  the  normal  distribution 
would  lead  to  acceptance.  Although  the  above  is  an  oversimpl i ri cat  ion , 
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it  does  indicate  the  possible  error  in  assuming  the  wrong  distribution 
in  determining  confidence  limits. 

The  third  source  of  error  in  assuming  constant  failure  rate  lies 
in  the  unrealistic  handling  of  suspension  data  in  estimating  population 
parameters.  Consider  as  an  examnle  problem  a  test  sample  of  failure 
times  given  as  3800,  3000,  4100  and  4200  and  a  sample  of  suspension 
times  given  as  3500,  4000,  4000  and  4500.  In  this  example,  eight  dif¬ 
ferent  components  were  tested  but  onlv  four  failed.  The  point  esti¬ 
mates  of  MTBF  assuming  exponential  and  Weibull  distributions  are  8000 
and  4130  respectively.  As  can  be  seen  the  'ITB*7  for  constant  failure 
rate  is  nearly  twice  the  value  for  the  nonconstant  failure  rate  Weibull 
distribution . 

In  the  first  two  sources  of  error  discussed  above  the  computed 
reliability  is  generally  conservative  for  components  and  systems  that 
wear  out.  Reliability  and  lower  confidence  limits  are  conservative  in 
that  lower  than  true  values  are  computed.  In  the  third  source  of  error 
the  resulting  reliability  is  generally  non conservative.  The  unrealis¬ 
tic  handling  of  suspension  data  seems  to  be  the  most  significant 
source  of  error  in  assuming  constant  failure  rate  for  testing  of 
mechanical  systems. 

A  second  assumption  or  requirement  which  is  often  made  in  relia¬ 
bility  analysis  is  that  the  data  sample  be  complete.  Consequently 
only  the  failure  times  in  a  given  test  are  considered  for  analysis. 

The  main  reason  for  this  assumption  is  that  statistical  methods  that 
treat  suspended  data  along  with  the  failure  data  are  in  many  cases 
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restrictive,  complicated  or  not  available.  It  is  clear,  however,  that 
suspended  data,  particularly  dnta  involving  suspension  times  that  are 
greater  than  some  or  all  of  the  failure  times  should  contribute  signifi¬ 
cant  information  toward  determining  the  unknovm  population  parameters. 
This  was  observed  in  the  above  discussion  of  the  constant  failure  rate 
assumption.  Suspended  data  can  result  in  manv  ways.  In  weapon  svst.em 
testing,  for  example,  it  can  arise  from  the  removal  of  a  component  from 
test  without  failure,  completion  of  a  system  tos^  prior  to  failure  or 
through  failure  resulting  directly  from  the  failure  another  compo¬ 
nent  or  cause  such  as  accident. 

A  third  assumption  which  is  often  made  in  determining  confidence 
limits  on  system  reliability  from  system  tests  is  that  the  system  can 
be  treated  as  if  it  were  a  single  component  with  no  differentiation 
being  made  as  to  which  component  within  the  system  actually  fails. 

All  failures  are  treated  as  svstem  failures.  A  necessary  underlying 
assumption  for  this  approach  is  that  either  all  components  have  con¬ 
stant  failure  rate  or  all  components  have  failed  a  number  of  times  so 
that  steady  state  conditions  prevail.  In  this  instance  confidence 
intervals  on  system  reliability  are  readi lv  derived  using  the  theory 
applicable  to  the  constant  failure  rate  case.  A  limitation  of  this 
assumption  is  that  results  of  individual  component  or  subsvstem  tests 
cannot  be  included  with  the  system  test  results.  Also,  if  components 
are  redesigned  during  the  course  of  a  test,  which  is  often  the  case  in 
large  weapon  system  testing,  total  test  time  on  the  redesigned  compo¬ 
nents  are  not  the  same  as  the  total  system  test  time.  The  usual 
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methods  of  analyzing  system  test  data  consequently  do  not  annlv  in 
this  case. 

Finally,  large  sample  theory  is  often  assumed  in  computing  confi¬ 
dence  intervals  since  methods  to  handle  small  samples  may  not  exist 
or  are  too  difficult  to  use. 

MAIN  FEATURES  AND  STATUS  OF  THE  PRESENT  MODEL 


The  reliability  data  analysis  model  presented  in  this  report  con¬ 
tains  a  number  of  general  solutions  to  overcome  many  of  the  limitations 
of  the  usual  assumptions  discussed  in  the  previous  section.  The  main 
features  of  the  model  are  summarized  as  follows: 

a.  Performs  goodness-of-fit  test  to  determine  the  best  fit 
probability  distribution  of  component  failure  times.  The  theoretical 
distributions  considered  are  the  exponential,  normal,  lognormal, 

Weibull  and  gamma. 

b.  Computes  maximum  likelihood  estimates  of  population  para¬ 
meters  for  general  theoretical  distribution  of  failure  times. 

c.  Can  handle  suspended  data  which  results  when  a  component  is 
tested  without  failure. 

d.  Computes  point  estimates  of  reliability  for  the  renewal  case; 
that  is,  for  the  non-constant  failure  rate  case. 

e.  Computes  lower  confidence  limits  on  component,  subsystem 
and  system  reliability  for  the  constant  failure  rate  case. 
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The  present  version  of  the  reliability  data  analysis  computer 


model  does  not  contain  all  of  the  features  ultimately  Planned  for  in 
the  final  version.  Most  significant  of  the  limitations  is  the  assump¬ 
tion  of  constant  failure  rate  in  determining  confidence  intervals 
although  point  estimates  reliability  are  determined  for  the  general 
non-constant  failure  rate  case.  In  addition  components  are  assumed  to 
be  in  series  for  system  reliability  computation  and  no  provision  is 
presently  made  for  preventive  maintenance  parts  replacement  of  compo 
nents.  The  analytical  methods  section  of  this  report  describes  the 
techniques  to  be  ennloved  in  removing  these  limitations  with  work  on 
future  versions  of  the  model  presently  being  undertaken. 

The  remainder  of  this  report  contains  a  discussion  of  the  analyti¬ 
cal  methods  used  for  computation  followed  hv  a  presentation  of  the 
computer  program  with  example  input  and  output  information. 


ANALYTICAL  METHODS 

The  purpose  of  this  section  is  to  briefly  outline  the  mathematical 
and  statistical  methods  used  or  planned  in  the  reliability  data 
analysis  model.  A  complete  discussion  of  probability  and  statistical 
theory  will  not  be  presented  in  this  report  and  the  reader  is  referred 
to  the  references  cited  for  more  complete  discussions.  Much  of  the 
theory  presented  is  straightforward  and  readily  found  in  the  litera¬ 
ture.  However,  some  of  the  statistical  and  computational  methods  used 
are  the  result  of  research  efforts  at  the  Watervliet  Arsenal  and  will 


IS 


be  the  subject  of  forthcoming  renorts  and  nubl i cations. 

There  are  a  number  of  good  texts  on  the  general  subiect  of  reli¬ 
ability.  The  texts  found  particularly  useful  to  this  writer  are  those 
by  Lloyd  and  Linow  [2],  Barlow  and  Proschan  [3],  Hnedenko,  Belyayev 
and  Solovyev  [4]  and  Pieruschka  [5]. 

Computation  of  Component  Mission  Reliability. 

Although  one  of  the  computational  goals  in  th^  model  is  system 
reliability  it  is  necessary  to  develop  the  analytical  methods  for 
computing  system  reliability  by  considering  the  individual  components 
or  elements  making  up  the  svstem.  This  is  particularly  true  in  the 
case  of  non-constant  failure  rate  components. 

At  the  present  tine  there  is  considerable  confusion  among  design 
engineers  on  the  definition  of  component  reliability  when  the  cormo- 
nent  is  part  of  a  system.  The  confusion  lies  primarily  on  the  time 
reference  used  in  computing  reliability.  Many  texts  on  reliability 
consider  time  in  terms  of  component  age.  However,  for  system  reliabi¬ 
lity  the  system  age  is  the  imnortant  time  reference.  Since  comnonents 
that  fail  within  a  system  are  replaced  or  renewed  at  random  points  in 

2 

Lloyd,  D,K.,  and  Linow,  M.,  "Reliability:  Management ,  Methods,  and 
Mathematics,"  Prentice-Mall,  Englewood  Cliffs,  New  Jersev  (1962). 

3 

Barlow,  R.E.,  and  Proschan,  F.,  "Mathematical  Theory  of  Reliability," 
John  Wiley  f,  Sons,  New  York  (1965). 

4 

Gnedenko,  B.V. ,  Belyayev,  Yu.K.,  and  Solovyev,  A.D.,  "Mathematical 
Methods  of  Reliability  Theory,"  Academic  Press,  New  York  (1969). 

’Pieruschka,  E.,  "Princinles  of  Reliability,"  Prentice-Hall,  Englewood 
Cliffs,  New  Jersey  (1963). 
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time,  component  apes  are  generally  not  Vnown  a  ^riori  as  a  function  of 
system  ape*  In  this  instance  renewal  theorv  must  ho  employed  fo 
determine  the  transient  mission  reliability  as  a  function  of  system 
age.  This  asnect  of  the  problem  will  he  considorod  in  the  following 
sections. 


(a)  Component  Reliability  Based  on  romnonent  Acq. 

Consider  first  a  single  con'^nen^  where  the  time  reference 
is  component  age.  Let 

f(t)  »  Probability  density  distribution  nr  coo^on^nt  failure 
times  or  tino  between  Failures, 
t 

F(t)  =  J  f(t)dt  =  cumulative  distribution  f motion. 

0 

=  Probability  that  the  coninnnent  will  fr  1  in  the  time 
interval  (0,t). 

R(t)  =  Reliability  defined  as  the  nrobahilitv  that  the  compo¬ 
nent  will  not  fail  in  the  time  interval  (0,t)  where 
the  component  is  assumed  new  at  time  0% 

=  1  -  F(t) 

Aft)  =  Conditional  failure  rate  whorft  A(t)dt  is  the  probability 
of  failure  in  the  time  interval  (t,t+dt)  "iven  that  the 
component  has  survived  to  time  t. 

The  conditional  failure  rate  >  ft/)  defined  above  is  a  useful  and 

descriptive  quantity  in  reliability  theory.  It  can  be  determined  from 

the  probability  distribution  of  failure  times  as  follows: 

Aft)  =  f(t)/fl-*(t)) 

=  f(t)/Pft)  fl) 

Figure  1  shows  descrintivelv  the  failure  ra*e  for  three  typical  cases 
in  reliability;  decreasing,  constant  and  increasing  failure  rate. 
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Component 

Failure 

Rate 


The  decreasing  failure  rate  implies  that  the  longer  a  component  survives 
the  lower  its  probability  of  failure.  This  is  typical  of  early  failures 
in  which  components  containing  manufacturing  or  material  flaws  tend  to 
fail  early.  The  increasing  failure  rate  is  typical  of  components  that 
wear  out  where  the  longer  a  component  survives  the  higher  the  probabi¬ 
lity  of  failure.  Mechanical  components  in  which  early  failures  have 
been  eliminated  generally  have  increasing  failure  rates  caused  by  such 
phenomena  as  fatigue,  corrosion,  erosion  and  abrasion.  The  constant 
failure  rate  is  typical  of  components  which  have  a  probability  of 
failure  independent  of  its  age.  This  can  result,  for  example,  from 
a  situation  in  which  excessive  loads  can  cause  failure  where  load 
fluctuations  are  purely  random  in  time  such  as  wind  loads. 
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Exarmies : 


1.  Exponential  Distribution 


For  this  case 


fft)  =  Ae 


-At 


-At 


F(t)  =  1-e 

R(t)  =  o"Xt 

A(t)  =  A  =  constant  f0') 

The  exponential  distribution  describes  the  constant  failure  rate  cas^. 
2 •  Weibul 1  Distribution 

ff  ,  -  t  l"1  e_ft/ri)B 


F(t)  =  1  -  e 


-ft/n) 


R(t)  =  e 


-ft/n) 


B  t 

Ut)  =  w  ^ 


B-l 


(V 


in  which  B  and  n  are  distribution  or  nonulation  naraneters.  Note  that 
for  the  Weibull  distribution  B  values  less  than,  eoual  to  and  greater 
than  unity  yield  failure  rates  which  are  decreasing,  constant  and 
increasing  resnectivelv.  This  characteristic  nakes  the  Weibull  distri¬ 
bution  a  useful  one  in  reliability  analysis. 


(b)  Mission  Reliability  of  a  Component  Within  a  System  [3, 

it  S]- 

Consider  next  the  situation  of  a  component  within  a  system* 
The  Quantity  of  interest  here  is  the  r^liabilitv  of  the  component  for 
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Barlow,  R.  E.,  and  Proschan,  F.,  "Mathematical  Theory  of  Reliability," 
John  Wiley  F,  Sons,  New  York  (1965). 
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Gnedenko,  B.  V.,  Belyayev,  Yu.K.,  and  Solovyev,  A.  D.,  "Mathematical 
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a  mission  time  interval  of  ft,  t  +  t)  where  t  is  the  system  age  and  t 
is  the  mission  length.  Prior  to  time  t  the  component  could  have 
failed  and  been  replaced  one  or  more  times,  figure  2  deoicts  the 
conditional  failure  rate  Aft)  of  a  component  within  a  svstem  showing 
failure  points  for  an  increasing  failure  rate. 


Component 
rai lure 
Rate 


Figure  2.  Component  Failure  Rate  vs  Svstem  Age  For  The  Case 

Of  Ideal  Repair 


In  general  the  failure  times  of  components  within  a  system  are  not 
known  in  advance  as  depicted  in  figure  2  and  consentient ly  must  be 
treated  probabilistically.  This  is  accomplished  by  defining  another 
quantity  hft)  called  the  renewal  rate  or  unconditional  failure  rate 
which  is  an  ensemble  average  of  the  failure  rate  over  the  population 
of  all  systems.  The  renewal  rate  is  defined  such  that 

h(t)dt  =  unconditional  probability  of  component  failure 
in  system  time  interval  (t,t+dt). 
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The  renewal  rate  for  a  given  component  is  a  function  of  the 
underlying  failure  distribution  given  by  the  following  equation: 


h(t)  =  f(t)  +  /r  h(x)f (t-x)dx. 

0 

Derivation  of  this  equation  is  nresented  in  references  [7],  [3]  ana 
elsewhere.  Logically,  equation  (4)  can  be  derived  from  the  theorems 
of  total  and  conditional  probabilities  [6]  along  with  the  definitions 
of  the  terms  in  equation  (4).  Multiplying  both  sides  of  eouation  (4) 
by  dt  and  considering  the  integral  as  a  sum  we  have 


f(t)dt 

=  Probability  that  the  original  component  fails 

in  the  time  interval  (t,  t  ♦  dt) . 

h (x)dx 

=  Probability  that  a  failure  and  subsequent  renewal 

occurred  in  time  interval  (x,  x  +  dx) . 

f (t-x)dt 

■  Conditional  probabilitv  that  a  component  which 

was  renewed  at  time  x  fails  in  the  time  interval 

(t,  t  +  dt). 

h  (x)f  ( t-x)dxdt 

=  Unconditional  probability  of  failure  in  time 

interval  (t,  t  ♦  dt)  for  a  component  which  could 

fail  at  time  x. 

f  (t)dt  + 


/^h (x) f ( t-x) dxdt 
n 

=  Total  probability  of  failure  which  is  the  sum  of 

all  possible  conditions  which  could  lead  to 

failure  in  the  time  interval  (t,  t  +  dt), 

=  h(t)dt  by  definition. 

The  mission  reliability  can  now  be  determined  from  the  renewal 
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rate  h(t)  using  the  relation  from  reference  [4] 


R(t,T) 

in  which 

R(t,T) 


l-F(t+x) 


1-F (t+T-x) 


h (x)dx 


1-F (t+T)  +  /t[l-F(t+T-x)]h(x)dx  (5) 

n 

Mission  reliability  at  svstem  time  t  for  a 
mission  length  t. 

Probability  of  no  failure  in  (t,  t+T) . 

Probability  that  the  original  component  in  the 
system  has  not  Failed  at  time  t+T. 

Probability  that  a  component  which  failed  at 
time  x  has  not  failed  at  time  t+T. 

Probability  of  failure  at  time  x. 


Typical  examples  of  the  renewal  rate  and  mission  reliability  of  an 
increasing  failure  rate  component  are  depicted  in  Figures  3  and  4.  As 
can  be  seen  the  renewal  rate  increases  with  svstem  age  until  about  a 
system  time  eaual  to  the  MTBf  (mean  time  betweep  failure)  of  the 
component.  The  renewal  rate  then  approaches  a  constant  value  enual  to 
the  reciprocal  of  the  ’1TBF.  The  reliability  on  the  other  hand  de¬ 
creases  from  a  value  of  1.0  at  t=0  and  approaches  a  constant  value. 

The  asymptotic  values  of  h(t)  and  R(t,r)  can  be  derived  directlv  from 
equations  (4)  and  (5)  by  passing  t  to  the  limit  infinity.  These 
values  are  given  by  the  relations 


h  (°°)  -  l/'ITBF 

1  00 

R(~,t)  =  U-F(x)]dx. 

T  (b) 


Gnedenko,  B.V.,  Belyayev,  Yu.K.,  and  Solovyev,  A.D.,  "Mathematical 
Methods  of  Reliability  Theory,"  Academic  Press,  New  York  (1069). 
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Component 

Renewal 

Rate 


mTbt 
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For  the  exponential  distribution  (constant  failure  rate)  the 

renewal  rate  and  reliability  are  both  constant  throughout  system  life 

-At 

and  are  equal  to  A  and  e  respectively. 

It  should  be  noted  at  this  point  that  for  a  high  reliability  sys¬ 
tem  (say  ,85  or  greater)  the  MTBF  of  the  components  comprising  the 
system  should  be  of  the  same  order  of  magnitude  as  the  expected  or 
required  system  life.  Consequently,  the  components  within  a  system 
are  essentially  exercised  within  the  transient  portion  of  their  lives 
thus  indicating  the  general  importance  of  considering  renewal  theorv 
in  determining  component  reliability  for  high  reliabilitv  systems. 

(c)  Numerical  Solution  for  the  Renewal  Rate. 

A  number  of  numerical  techniques  were  investigated  to 
solve  the  integral  equation  (4)  for  the  renewal  rate  h(t).  The  most 
computationally  efficient  solution  investigated  thus  far  was  through 
the  use  of  finite  difference  methods  [7] .  finite  difference  yields  a 
direct  solution  for  the  renewal  rate.  Other  more  general  technioues 
investigated  which  give  complete  solutions  to  the  renewal  problem 
(e.g.  probability  distribution  of  total  number  of  failures  or  time  to 
nth  failure)  involved  the  use  of  orthogonal  expansions  of  Laplace 
transforms.  Two  sets  of  orthogonal  functions  considered  were  trigono¬ 
metric  and  Laguerre  polynomials.  The  finite  difference  solution, 
although  limited  to  solution  of  the  renewal  rate,  generally  reouired 
less  computer  time  for  given  accuracy.  The  solution  is  also  quite 
general  in  that  it  can  be  used  to  determine  renewal  rate  for  most 
theoretical  distributions  of  failure  time  applicable  to  reliability 

7 

McKelvey,  R.  W.,  "An  Analysis  of  Approximate  Methods  for  Fredholm’s 
Integral  Equation  of  the  First  Kind,”  December  1056,  AD650530. 
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theory. 

In  solving  equation  (4)  using  finite  differences,  the  density 
f(t)  and  the  renewal  rate  h(t)  are  discretized  over  a  fixed  time  inter¬ 
val  generally  the  system  life.  In  this  case  f(t)  and  h(t)  are  written 
as 


f(t) 

=  f (kAt)  =  iy 

h  (t) 

*  h  (kAt)  =  h^ 

t 

=  kAt  where  k=0,l 

The  integral  equation  (4)  can  then  he  written  as 


Jf(t-x)h(x)dx  «  —  l  [f (kAt-iAx)h(jAx)  ♦ 

0  2  i„o 

♦f  (kAt-  (j+l)Ax)h  ( (i  +  l)Ax)  ] 


Ax 


k-1 


ox  y 
2  L 


j  =  l  [f,  .h.  ♦  f,  .  h.  .] 

J  k-i  j  k-i -1  i  +  lJ 


for  k  =  1 , _ ,n 

0  for  k  =  o. 


(7) 


Substituting  into  equation  (4)  then  gives 


h 


k 


Ax 


k-1 


fk  +  2  l  ffk-ihi  *  fk-i-lVll 

j=0  •  ' 


Tor  k  =  1 , . . .  ,n 
for  k  =  0. 


(S') 


Equation  (8)  renresents  n+1  equations  with  n+1  unknowns  which  can  be 
readily  solved  using  Gaussian  elimination. 
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A  difficulty  is  encountered  in  some  situations  where  f(0)  =  00 
such  as  in  the  case  of  the  Weibull  distribution  with  shape  parameter 
8<1.  In  this  instance  f(0)  is  fixed  at  a  finite  value  such  that  the 
actual  area  under  f(t)  in  the  interval  (0,  At)  (i.e.  F(At))  is  made 
equal  to  the  finite  difference  area  (f(0)  +  f(At))At/2. 

(d)  Numerical  Solution  for  Mission  Reliability. 

Once  the  renewal  rate  is  determined  the  transient 
mission  reliability  can  then  be  computed  from  equation  (5): 

R(t,x)  =  l-F(t+x)  +  /t  fl-^(t+T-x)]h(x)dx. 

n 

Numerically,  this  equation  is  difficult  to  solve  since  the  integral 
must  be  evaluated  over  the  whole  interval  (0,t)  for  each  different 
value  of  t.  By  making  use  of  the  renewal  equation  an  eouivalent 
equation  for  R(t,x)  can  be  derived  which  is  computationally  more 
efficient.  Integrating  the  renewal  equation  (4)  from  0  to  T  gives 

.T 

H(T)  «  F (T)  +  /  H(T-x)f(x)dx  (9) 

0  j 

where  H(T)  =  /  h(t)dt. 

0 

Making  the  substitution  T=t+x  in  equation  (9)  gives  the  following 
equation : 

t+x 

H(t+x)  =  F (t+x)  +  /  H(t+x-x)f (x)dx.  (10) 

0 

Integrating  by  parts  and  using  the  fact  that  H(0)  =  0  and  FfO)  =  0 
then  gives 


H(t+T) 


t+T 

=  F  (t+T )  +  /  P  fx)h  (t+T-x)dx 
0 

t+T 

=  r(t+T)  +  /  17  (t+T-x)h  (x)dx.  (11) 

0 

From  equation  (11)  the  following  eouation  can  he  readilv  derived: 
t  t+T 

/  [l-F(t+T-x)]h(x)dx  =  r (t+x)  -  J  [l-r(t+T-x)]h (x)dx 
0  t 

(12) 

t 

where  the  definition  H(t)  =  /  h(x)dx  was  used* 

0 

Substituting  equation  (12)  into  equation  (5)  for  mission  reliability 
finally  yields 

t+x 

R(t,T)  =  1  -  /  [1-F (t+T-x) ]h (x)dx  (13) 

t 

Note  that  in  this  equation  the  integral  is  evaluated  onlv  over  the 
interval  (t,  t+T)  and  that  the  value  of  F(y)  is  renuired  only  over 
the  interval  (0,t).  This  simplifies  the  numerical  solution  for  mission 
reliability. 


( e )  Average  Mission  Reliability. 

Interval  reliability  as  discussed  in  the  previous 
section  is  a  transient  function  of  svstem  age.  Weapon  svstem  reouire- 
ments,  however,  generally  specify  one  value  o^  mission  reliability  for 
the  system.  This  specified  value  could  represent  the  lowest  reliability 
to  be  experienced  by  the  system  or  an  average  value  over  system  life* 

In  the  computer  model  average  system  reliability  is  computed  using  the 


relation 


(14] 


in  which 

Rj (r)  =  Reliability  for  the  ith  mission 

n  =  Exnected  number  of  missions  over  system  life. 


Comnutation  of  System  Mission  Reliability. 


In  general,  system  reliability  can  be  computed  directly  from 


the  component  reliabilities  using  an  appropriate  reliability  model 
[2,8],  For  example,  for  the  series  reliability  model,  system  reliabili¬ 
ty  Rs (t)  is  determined  from 


n 


(15) 


i=l 


in  which 

R^(t)  3  Reliability  of  the  ith  component 

n  =  Total  number  of  components. 

A  series  model  is  applicable  whenever  the  failure  of  a  single 
component  within  a  system  results  in  a  failure  of  the  system.  This 
model  has  been  initially  chosen  for  the  data  analvsis  computer  model. 

In  the  more  general  case  where  redundancy  and  load  sharing  are 
inherent  in  the  system,  reliability  models  are  more  complicated  but 
can  be  derived  in  most  given  situations  [2,8]. 
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Bazovskv,  I.,  "Reliability  Theorv  and  Practice,"  Prentico-IIal 1 , 
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Probability  Distribution  Selection. 


One  of  the  more  difficult  and  questionable  aspects  of  data 
analysis  is  the  selection  of  a  theoretical  distribution  applicable  to 
a  given  set  of  data.  This  is  particularly  true  for  small  data  samples 
where  information  is  ultimately  renuired  at  the  tails  of  the  distri¬ 
bution.  In  the  case  of  component  failure  data,  theoretical  considera¬ 
tions,  prior  history  and  experience  play  a  large  part  in  distribution 
selection.  For  example,  for  particular  failure  modes  such  as  fatigue 
the  lognormal  and  Weibull  distributions  have  been  found  to  yield  good 
characterizations  of  the  data  [9,  10], 

The  data  analysis  computer  model  includes  a  distribution  fitting 
program  to  assist  in  the  selection  of  a  best-fit  theoretical  distri¬ 
bution  for  use  in  reliability  computation.  Essentially,  this  program 
computes  the  standard  error  for  a  number  of  different  candidate 
distributions.  The  standard  error  is  a  measure  of  the  deviation  of 
the  data  from  the  theoretical  cumulative  distributions. 

The  Kolmogorov-Smimov  statistical  goodness-of-fit  test  is  also 
used  to  determine  which  of  the  theoretical  distributions  can  be 
rejected  for  given  confidence  level.  Final  selection  of  the  distri¬ 
bution  can  then  be  made  based  on  the  computer  results,  theoretical 
considerations  and/or  personal  experience  and  judgment. 

(a)  Candidate  Distributions. 

Following  are  the  candidate  distributions  considered 
for  characterizing  component  failure  times.  In  order  to  standardize 
terminology  used  to  describe  the  distribution  parameters,  the  terms 
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^Dolan,  T.  J.,  "Basic  Concents  of  Fatieu*  Damage  in  Metals,"  Metal 
Fatigue,  McGraw-Hill,  New  York  (1959),  Chanter  3,  on.  3  -67. 
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Freudenthal,  A.  M. ,  and  Shinozuka,  M.,  "Structural  Safety  Under 
Conditions  of  Ultimate  Load  Failure  and  ^atipue,"  Columbia  Universitv 
New  York,  N.Y.,  October  1961,  WADD  Technical  Reoort  6] -177, 


scale,  shape  and  location  parameters  are  used  for  all  distributions 
and  are  defined  below. 


where 


where 


where 


(1)  Exponential 


f(t)  = 

l  -t/n 

—  e 

n 

n  = 

MTBF  =  scale  parameter 

=5 

l/X 

X  = 

constant  failure  rate. 

Normal 

1  t-P  2 

f(t)  - 

_L  -7(— > 

e 

/2u  a 

P  = 

KfTBP  =  scale  narameter 

a  a 

Standard  deviation  =  shar>e 

Weibull 

f(t)  = 

3  t-Y  e-1  -(^)6 

(p)  ("7p  e 

P  = 

2 

a  = 

‘itbf  =  y  +nr(i+i/B) 

*>  2 

Variance  =  rf  [r(l+2/B)-r  (1 

n  = 

Scale  parameter 

3  = 

Shane  parameter 

Y  * 

Location  parameter. 

Lognormal 

f(t)  » 

1  -  ^^(£n  (t-Y) -T)) 

(t-Y)^2rF  3 

(16) 


(17) 


(18) 


(19) 
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where 

\i 

=  >rr  bf  = 

w(n*B2/2) 

2 

a 

=  Variance 

.  (/-l) 

n 

=  Average 

(in  t]  =  Scale  parameter 

B 

=  Standard 

Deviation  (&n  t)  =  Shape  parameter 

Y 

=  Location 

parameter. 

(5)  Hamma 

f(t)  = 

8-1 

(t-y) 

r(8)n6 

U  = 

*1TBF  = 

y+8h 

a2  = 

Variance 

-  Bn2 

n  = 

Scale  parameter 

B  = 

Shape  parameter 

Y  - 

Location 

parameter 

More  detailed  discussions  of  these  distributions  can  be  found 
in  the  literature  [6,11], 


(b)  Least  Squares  Fit  of  Data  [12]. 

For  each  distribution  other  than  Samma  a  least 
squares  fit  of  the  data  to  the  distribution  is  made.  Henerallv,  the 
theoretical  distribution  is  linearized  as  far  as  possible  to  simplify 
computation.  In  the  case  o*  the  Samma  distribution  maximum  likelihood 
estimates  of  parameters  are  used  rather  than  solving  the  more  difficult 
nonlinear  problem. 

The  cumulative  distribution  function  F(t)  is  used  in  fitting  the 
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data  with  the  median  ranks  being  used  as  the  true  values  of  the  cumula¬ 
tive  distribution  associated  with  each  ordered  data  noint  [13],  The 
median  ranks  for  a  complete  data  sample  are  computed  using  the  relation 


j-0.3 


n-0.4 


(21) 


where 


n  =  Data  sample  size 

i  -  1 , . . .  ,n 

=  Order  numbers  for  the  data  where  the  data  is 
sorted  in  increasing  order. 


When  suspended  data  has  been  generated  along  with  failure  data, 
the  median  ranks  are  determined  for  the  failure  data  onlv  with  the 
suspension  times  being  used  to  modifv  the  median  ranks  fl  3] .  In  this 
instance  the  failure  and  suspension  data  are  sorted  in  increasing 
order.  The  sample  size  n  in  equation  (21)  is  now  the  total  number  of 
failure  and  suspension  data  items.  The  order  number  i  associated  with 
each  failure  item  is  determined  as  follows;  initiallv  i  is  set  equal 
to  zero.  An  increment  is  then  computed  which  is  to  be  added  to  i  to 
give  the  order  number  for  the  first  failure  item  using  the  general 
relation 


New  Increment 


n  ♦  1  -  (Previous  Order  Number) 


1+ (Number  of  Items  nol lowing  the  Present 
Suspension  Set). 


(22) 


If  there  are  no  suspension  data  prior  to  the  first  failure,  then  the 
increment  is  1.0  and  the  order  number  for  the  first  failure  item  is 
1,0.  This  procedure  is  repeated  for  each  subsenuent  failure  item 


13 


'Johnson,  L.O.,  "Theory  nnd  Technique  of  Variation  Pesenrch,"  HI  sevier, 
New  York  (1964),  Chapter  8.  ^ 


using  equation  (22)  in  each  case.  The  median  ranks  are  then  computed 
using  equation  (21). 


The  median  ranks  are  considered  to  be  the  true  values  of  ^(t) 


associated  with  the  ordered  failure  times  t^.  The  distribution  para¬ 
meters  for  a  given  theoretical  form  of  F(t)  are  then  determined  by 
minimizing  the  sum  of  the  squares  of  the  error  between  the  theoretical 
distribution  evaluated  at  the  failure  times  and  the  median  ranks  or 
between  the  transformed  distribution  and  median  ranks. 

Consider  as  an  example  the  Weibull  distribution.  In  this  case 

t-Y  8 

F(t)  =  1  -  e  .  (23) 

Rearranging  equation  (23)  and  then  taking  double  logs  gives 


inin(l/(l-F(t)))  *  B (in (t-y) -Ann) .  (24) 

For  fixed  Y,  this  equation  is  linear  in  £n(t-Y).  The  error  between 
the  transformed  theoretical  distribution  and  the  data  is  then  given  hv 

ei  *  B(An(t.-Y)  -  inn)  -  inin(l/(l-MR.))  (23) 

where  t.  =  ith  failure  time 

l 

MR.  =  ith  median  rank. 

l 

The  parameters  B,  n  and  y  are  then  found  which  minimize  the  total 
2  2 

square  error  e  =  A  similar  approach  is  used  for  other  distri¬ 

butions  . 


(c)  Computation  of  the  Standard  Error  [12], 

In  general,  the  standard  error  in  the  fitting  of 


^Draper,  N.R.,  and  Smith,  H,,  "Applied  Regression  Analysis"  John  Wiley 
§  Sons,  New  York  (1966). 
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a  theoretical  distribution  to  data  is  not  the  same  as  the  error  dis¬ 
cussed  in  the  previous  section.  The  standard  error  is  a  measure  of 
the  difference  between  the  untransformed  theoretical  distribution  and 
the  median  ranks.  The  standard  error  is  defined  bv  the  following 
equation: 

(  l  (F(t.)-MR  )2)1/2 

i  =  l  1  1 

s.e.  =  Standard  Error  =  - —  (76) 

n  -  n  -  1 


where 


P(t.) 


Theoretical  distribution  evaluated  at  the  failure 
tines  t^. 


MR^  »  Median  ranks 

n  =  Number  o^  failure  points 

p  =  Number  of  population  parameters  estimated  in 
determining  R(t). 

The  smaller  the  s.e.  the  closer  the  data  fits  the  theoretical  distri¬ 
bution. 


(d)  Kolmogorov-Smimov  Goodness  of  Fit  N4]. 

A  non-parametric  distribution  has  been  derived  bv 
Kolmogorov  f 1 5 ]  and  Smirnov  f 1 6]  for  a  particular  statistic  d  which 
is  a  measure  of  the  fluctuation  of  sample  data  about  the  theoretical 
distribution  from  which  the  sample  is  drawn.  The  statistic  d  is  de¬ 
fined  as  the  maximum  absolute  deviation  between  the  theoretical  and 
observed  cumulative  probability  distributions.  The  theoretical  distri¬ 
bution  is  fixed  by  specifying  both  the  functional  form  and  the  para¬ 
meter^ _ 

A*Siegel,  S.,  "Nonparametric  Statistics  for  the  Behavioral  Sciences,” 
McGraw-Hill,  New  York  (1956),  Chapter  4. 


39 


5Kolmogorov,  A.,  "Confidence  Limits  for  an  Unknown  Distribution 
Function,"  Annals  Mathematical  Statistics,  Vol.  12,  (1D41),  nn.  461- 
463. 
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In  application  of  the  K-S  statistic,  a  given  set  of  data  is 
hypothesized  to  have  been  drawn  from  a  given  theoretical  distribution. 
The  statistic  d  is  then  determined  from  the  data  and  compared  to  the 
K-S  distribution  of  d.  That  is,  for  given  significance  level  a  or 
confidence  level  (1-a)  the  theoretical  value  d^  is  determined  from 
K-S  tables  and  compared  to  the  computed  d.  If  d>d^  ,  the  hypothesis 
is  rejected.  If  d<d^  then  there  is  no  sufficient  reason  to  reject  the 
given  theoretical  distribution  and  the  hypothesis  may  be  accepted. 

The  statistic  d  is  determined  for  the  candidate  distributions 
in  the  computer  model  to  provide  bases  for  rejecting  distributions  and 
to  help  indicate  the  best-fit  distribution. 

Maximum  Likelihood  Estimates  of  Parameters  fl7]. 

Maximum  likelihood  estimates  of  population  parameters  are 
determined  for  each  component  within  a  svstem  for  which  failure  data 
has  been  generated.  There  are  a  number  of  reasons  for  generating 
maximum  likelihood  estimators  as  part  of  the  reliability  model.  First, 
maximum  likelihood  yields  estimates  of  parameters  which  have  a  number 
of  desirable  attributes.  For  example,  if  an  efficient  estimator  for 
small  samples  exists  (i.e.  one  with  minimum  variance),  then  maximum 
likelihood  provides  such  an  estimator.  In  general  then  minimum  con¬ 
fidence  intervals  can  be  derived  in  this  case.  Maximum  likelihood 
estimators  are  also  consistent  in  that  thev  approach  the  true  para¬ 
meter  values  as  sample  size  increases.  In  addition,  maximum  likeli¬ 
hood  estimators  are  asymptotically  normal  as  sample  size  increases 

17 

Mood,  A.M. ,  and  Craybill,  V.A.,  "Introduction  to  the  Theory  of 
Statistics,"  McGraw-Hill,  New  York  (1963),  Chapter  8. 
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for  most  cases.  This  simplifies  determination  of  confidence  intervals 
for  large  samples. 

A  second  reason  for  determining  maximum  likelihood  estimates  is 
that  they  can  be  computed  for  relatively  general  underlying  probability 
distributions  through  the  use  of  routine  mathematical  procedures. 
Finally,  the  numerical  procedures  planned  to  be  used  in  the  computer 
model  to  determine  confidence  intervals  require  the  maximum  likelihood 
estimates  of  parameters  to  improve  computational  efficiency. 

The  basic  idea  behind  maximum  likelihood  estimation  is  relatively 
straightforward.  One  assumes  first  that  a  sample  of  n^.  failures  and 
ns  suspensions  have  been  generated  from  a  given  theoretical  distribu¬ 
tion  with  parameters  *  (ct^  . . .  .ot^)  .  The  failure  and  suspension 
data  are  designated  as  t^  (i=l , . . .  ,n£.)  and  t  .  (i=l , . . .  ,n^)  respectively. 
The  joint  probability  distribution  of  the  random  sample  of  failure  and 
suspension  times  can  be  written  as 


L(t_;a) 

3 

gt(tfl,..*,tfh  ,tsl,...,tsn  ;a) 

—  Is 

where 

L(t;a) 

- 

Defined  as  the  likelihood  function 

8i 

- 

Joint  distribution  of  the  sample  outcome 

t_ 

3 

Vector  of  the  sample  data  values 

3 

Vector  of  parameters  for  given  underlying 
distributions  of  failure  tines. 

It  is  next  assumed  that  each  failure  and  suspension  time  is  a 
statistically  independent  outcome.  Equation (27)  can  then  be  written 
in  terms  of  the  failure  density  f(t)  and  the  cumulative  distribution 
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F (t)  as  follows  [6] 


nf 


L(l;o)  =  C  |  |  f(tfi;«)  |  |  fl-F(t  -a)]  (28) 

i=*l  i  =  l  J 


where  C 

f (t ;a) 
F (t ;a) 


Normalizing  constant  such  that  the  area  under 
L(tja)  is  unity 

Theoretical  density  distribution  of  failure  times. 
Cumulative  distribution  function 


=  Jf  f(x;nt)dx 
0 

[1-F(t;a)]  -  Probability  that  the  failure  time  is  greater  than 

t.  This  term  is  used  to  represent  the  nrobabilitv 
of  obtaining  the  suspension  times. 


As  can  be  seen,  the  likelihood  function  L  is  defined  as  the 
multivariate  probability  distribution  of  the  random  failure  and  sus¬ 
pension  times.  Solving  for  the  parameters  such  that  L  is  maximized 
consequently  gives  the  highest  probability  density  for  the  given  sample 

A 

outcome  t.  The  parameter  values  which  maximize  L  are  denoted  as  a 
and  are  called  the  maximum  likelihood  estimates  of  ou 

In  practice  Jin  L  given  by  equation  (29)  is  maximized  rather  than 
L  to  simplify  computation.  This  can  be  done  since  the  maximum  of 
any  positive  function  and  its  log  are  equivalent. 


Jin  L(t_;a) 


Jin  C  ♦ 


in  f(tfi;ct) 


♦ 


l  in  f 1-F (t  .;a)l 
j-1  “ 


(29) 


^Panoulis,  A.,  "Probability,  Random  Variables ,  and  Stochastic 
Processes,"  McGraw-Hill,  New  York  (196S). 
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Consider  the  Weibull  distribution  as  an  example.  For  this  case 


6-i  r— i 

t-Y*  "I— J 1 


f(t;a)  «  -(— ) 

F(t;ct)  =  1  -  e  n 


where 


a  ■  (n,B,Y). 


The  likelihood  function  can  then  be  written  as 

nf 

&n  L  *  £n  C  +  ][  [£nB-£nn+ (B-l) (£n(t_.-y) 

i=l  fl 

"s  tsi-YS 

-  *nn)  -  (-— )  ]  -  l  (—£—■)  (30) 

n  j-i  n 

Solving  for  n , B  and  y  which  maximize  % n  L  vields  the  maximum  likeli¬ 
hood  estimators  for  the  Weibull  parameters. 

Likelihood  functions  such  as  given  by  equation  (30)  are  maximized 
in  the  computer  model  using  Rosenbrock's  algorithm  f 1 8]  which  is  a 
general  solution  of  the  unconstrained  minimization  problem  for  non¬ 
linear  functions.  Minimizing  the  negative  of  £n  L  maximizes  the 
function. 

Once  point  estimates  of  the  failure  distribution  parameters  are 
generated  for  a  given  component,  the  mission  reliability  can  be  com¬ 
puted  using  the  methods  previously  presented. 

18Rosenbrock,  H.H.,  "Automatic  Method  for  Finding  the  Greatest  or 
Least  Value  of  a  Function,"  Computer  Journal,  Vol.  3,  (1^60), 
pp.  175-184. 
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Bayesian  and  Classical  Confidence  Intervals  for  Component 


Reliability. 

Classical  methods  of  confidencing  component  mission  reli¬ 
ability  for  small  samples  are  generally  available  onlv  for  the  constant 
failure  rate  case.  In  this  instance  the  y  distribution  reoresents 
the  inferencing  distribution  for  the  ‘fTB^  from  which  confidenced  reli¬ 
ability  can  be  derived  [2],  Bayesian  methods  will  consequently  be  used 
the  computer  model  for  the  non-constant  failure  rate  comnonents  since 
this  method  provides  a  systematic  means  of  determining  confidence  inter¬ 
vals  for  more  general  problems.  A  complete  discussion  of  Bayesian 
statistics  will  not  be  presented  in  this  report  and  the  reader  is 
referred  to  the  books  by  Lindley  [19]  for  a  more  complete  discussion 
of  this  topic. 

Bayesian  inferencing  is  based,  as  the  name  implies,  on  Bayes' 
probability  theorem  which  is  essentially  a  conditional  probability 
statement : 


Bayes*  Theorem: 


f(y,x)  =  f(v;x)f(x)  =  f(x;y)f(y) 

in  which 


[51) 


f(y,x) 

f (y;x) 
f(x;y) 


Joint  probability  density  of  random  variables 
x  and  y. 

Conditional  probability  of  y  given  x 
Conditional  probability  of  x  given  y 


Lloyd,  D,  K.,  and  Lipow,  M. ,  "Reliability:  ‘Management,  Methods,  and 
Mathematics,"  Prentice-Hall,  Englewood  Cliffs,  New  Jersey  (1962). 
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Lindley,  D.V.,  "Introduction  to  Probability  and  Statistics  from  a 
Bayesian  Viewpoint,"  Part  1:  Probability  and  Part  2:  Inference, 
Cambridge  University  Press,  Cambridpe  (1965). 
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fOO,  f(y)  = 

Marginal  distribution  of  x  and  v  respective! v. 

From  equation  (31) 


f(y;x)  = 

f(y) 

f(x;y)  t—  *  OP 

f(x) 

In  the  Bayesian 

approach  the  unknown  Parameters  of  the  distribu- 

tion  of  failure  times  are  considered  to  he  random  variables.  The 
sample  outcome  t  from  a  given  test  arc  also  random  variable  which 
are  dependent  on  the  distribution  Parameters %  Bayes’  theorem,  equation 


(32),  can  be  written 

in  this  instance  as 

f(a;x)  = 

f(a) 

f(x;oO  -  03) 

f(x) 

in  which 


f(a;x)  = 

Density  distribution  of  parameters  given  the 
test  sample  x. 

f(x;a)  = 

Density  distribution  of  the  sample  outcome 
given  the  parameters  ot. 

Marginal  densities  of  nt  and  x  respectively. 

In  equation  (33)  the  distribution  f(x;nt)  is  lust  the  livelihood 
function  L(xja)  hy  definition.  The  distribution  f(nt)  is  called  the 
prior  of  a  and  f(cx;x)  is  called  the  posteriori  distribution  of  n 
since  it  is  determined  after  the  sample  is  determined  and  is  con¬ 
sequently  conditioned  by  the  test  results. 

The  prior  distribution  f(a)  generally  reflects  prior  knowledge 
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or  prior  data  on  the  nooulation  parameters  a.  before  the  test  results 
x  have  been  generated.  For  the  case  of  weapon  system  components  there 
is  at  the  present  time  little  or  no  prior  information  and  consequently 
f(ct)  must  reflect  in  some  way  our  orior  ignorance  of  the  values  of  a. 
Choosing  the  prior  distribution  to  reflect  degrees  of  ignorance  is  one 
of  the  more  difficult  and  questionable  aspects  of  Bayesian  statistics. 
Much  research  work  is  presently  being  conducted  on  the  question  of 
prior  distribution  selection.  Priors  which  represent  maximum  ignorance 
for  some  of  the  simpler  problems  have  been  found  [10,70] .  It  has  also 
been  found,  however,  that  a  uniform  prior  on  the  parameters  yields 
confidence  intervals  that  are  generally  close  to  exact  in  a  classical 
frequency  sense  for  a  number  of  reliability  indices.  The  uniform 
distribution  assigns  equal  probability  to  all  values  of  a  random 
variable  over  a  given  range.  The  robustness  of  uniform  priors  is 
discussed  somewhat  in  reference  [19].  Uniform  priors  have  consequently 
been  used  in  our  reliability  statistics  work  and  have  also  been  found 
to  be  generally  robust. 

The  distribution  f(x)  in  eauation  (33)  is  considered  to  be  a 
constant  and  is  determined  so  that  the  area  under  ffojx)  i-s  unity. 

Equation  (33)  can  now  be  rewritten  in  the  following  form: 

f(a;x)  »  C  L(xja)  (34) 

where  f(a;x)  =  Posteriori  distribution  of  given  x_. 

C  =  Constant  such  that  area  under  f(3>x)  is  unity. 

The  constant  C  contains  the  constant  terms 
f(x)  and  f (a) . 


19 

Lindley,  D.  V.,  "Introduction  to  Probability  and  Statistics  crom  a 
Bayesian  Viewpoint,"  Part  1:  Probability  and  Part  2:  Inference, 
Cambridge  University  Press,  Cambridge  (1905). 

20 

Jaynes,  E.  T.  ,  "Prior  Probabilities,"  IEEE  Transactions  on  Svste^s 
Science  and  Cybernetics,  Vol.  SSC-4,  No.  3,  Sent.  1^68,  nn.  227-241. 
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L(x^;a)  =  Likelihood  function. 

In  computing  confidenced  reliability  our  interest  is  not  on  the 
parameters  themselves  but  rather  on  a  function  of  the  parameters. 

From  probability  theory  [6]  the  cumulative  distribution  of  any  function 
z  of  the  random  variables  a  is  determined  by  the  relation 


Fz(z;x) 

in  which 

Vz(z;x) 

Dz 

f(a;x) 


///  f(a»x)<ta 


(35) 


Posteriori  distribution  of  Z 

Domain  of  I_  such  that  ZQx)<z 

Posteriori  distribution  of  ft  given  by 
equation  (34). 


Once  the  posteriori  distribution  n(z;x)  is  known  confidence 
intervals  can  be  constructed  on  z  for  any  Riven  confidence  level  CL 
by  determining  the  annronriate  Percentage  points  directlv  from  F(z;x^. 
For  example,  a  lower  confidence  limit  is  determined  by  solving  the 


following 

equation 

for 

ZL: 

FZ ^ZL’— ^ 

se 

1  -  CL 

in  which 

ZL 

3 

Lower  confidence  limit. 

CL 

a 

Confidence  level. 

Since  component  mission  reliability  is  a  direct  function  of  the 
failure  distribution  and  hence  of  the  parameters,  confidence  intervals 

Panoulis,  A.,  "Probability,  Random  Variables,  and  Stochastic 
Processes,"  McCraw-Hill,  New  York  (1965). 


50 


on  mission  reliability  can  be  determined  using  the  above  formulation* 
Solution  of  equations  (34),  (35)  and  (36)  for  confidence  limits  can 
be  accomplished  numerically  using  comnuter  routines. 

Bayesian  Confidence  Intervals  for  System  Reliability. 

(a)  Constant  failure  Rate  Components  In  Series. 

Confidence  intervals  on  component  reliability  for  con- 
stant  failure  rate  can  be  determined  using  the  X"  distribution  [2], 

For  a  system  of  components  in  series,  however,  classical  confidence 
intervals  can  be  readilv  determined  only  if  all  components  have  enual 
test  times.  In  this  instance  the  system  can  be  treated  as  if  it  were 

a  component  with  all  failures  being  counted  as  system  failures  regard- 

2 

less  of  which  components  fail.  The  y  distribution  can  again  be  used 
to  determine  confidence  limits. 

For  series  components  with  unequal  test  times,  however,  classical 
methods  do  not  generally  apply.  For  this  case  Bayesian  statistics,  as 
discussed  in  the  previous  section,  can  be  employed  to  determine  confi¬ 
dence  intervals. 

Consider  the  series  system 


R  =  M  R-  (37) 

S  .1 

i=l 

in  which  R  =  Svstem  reliability, 

s 

n£  =  Number  of  components. 

R.  =  Reliability  of  the  ith  component, 
l 

2 

Lloyd,  D.  K, ,  and  Linow,  M. ,  "Reliability:  Management,  Methods,  and 
Mathematics,"  Prentice-Hall,  Englewood  Cliffs,  New  Jersey  (1962). 
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For  constant  failure  rate,  =  e  ^iT  and  the  system  reliability 
becomes 


R 

s 


i  =  l 


-X  T 
s 
e 


where 

X 

s 


n 


c 

l  K 


i=l 


(38) 


The  system  failure  rate  X  is  eou<al  to  the  sum  of  the  individual 

s 

component  failure  rates.  The  posteriori  probability  density  of  the 

individual  component  failure  rates  can  be  determined  using  equation 

(34)  by  letting  n  =  X  .  Since  X  is  equal  to  the  sum  of  independent 
—  i  s 

random  variables,  the  distribution  for  X  is  equal  to  the  convolution 

s 

of  the  individual  component  distributions  [6].  Performing  the  required 
convolution  is  generally  tedious  and  it  is  at  this  point  that  a  simpli¬ 
fying  assumption  is  made.  Using  the  Central  Limit  Theorem  for  sums  of 
random  variables  [6]  it  is  assumed  that  the  posteriori  distribution  for 
X$  is  Gaussian,  Consequently,  all  that  is  required  is  the  mean  and 

variance  of  X  to  define  the  posteriori  distribution.  These  are  deter- 
s 

mined  using  the  relations 


E(Xs) 

Var(A$) 

Where  E(*) 

Var(*) 


"c 

l  E(X.) 
i=l 

nc 

l  Var (X. ) 
i=l 

Expected  value 
Variance 


(30) 


(40) 


^Papoulis,  A.,  "Probability,  Random  Variables,  and  Stochastic 
Processes,"  McGraw-Hill,  New  York  f1Q65). 


Confidence  intervals  can  now  be  constructed  for  Xs  and  hence  for  system 
reliability  using  equations  (39)  and  (40). 

The  accuracy  of  making  the  Gaussian  assumption  was  checked  and  it 
was  found  that  a  system  with  about  ten  components  and  one  to  three  fail¬ 
ures  experienced  ner  component  yielded  relatively  accurate  confidence 
intervals  in  comparison  to  true  values.  In  the  example  problem  consi¬ 
dered  test  times  for  each  component  were  assumed  equal  which  permitted 
the  computation  of  classical  confidence  limits. 


(b)  Non-constant  Failure  Rate  Components  In  Series, 
For  the  series  svstem  conation  (37)  applies: 


R  =  |  |  R 

s  i  =  l 

TakinR  logs  of  both  sides  of  this  enuation  yields 


In  R  =  )'  InR. 


i=l 


Equation  (41)  represents  a  sum  of  independent  random  variables 
and  consequently  the  posteriori  distribution  of  lnR^  can  be  determined 
using  the  Gaussian  assumption  as  discussed  in  the  previous  section. 
Confidence  intervals  can  then  be  constructed  from  this  distribution. 


COMPUTER  PROGRAM 

Two  computer  programs  have  been  prepared  to  perform  the  compu¬ 
tations  required  for  reliabilitv  data  analysis.  One  program  called 
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DISTSEL  performs  goodness-of-fit  computations  and  is  used  to  assist  in 
distribution  selection  given  component  failure  and  suspension  data. 
Included  in  this  program  are  computations  of  the  standard  error  and 
Kolmogorov-Smimov  statistic.  At  the  present  time,  the  theoretical 
distributions  that  can  be  handled  are  the  exponential ,  two  and  three 
parameter  Weibull,  two  parameter  lognormal,  two  parameter  namma  and 
the  normal . 

The  second  commuter  program  called  RELIAB  performs  commutations 
required  to  determine  reliability  and  confidence  limits  for  components, 
subsystems  and  system  given  component  data.  The  component  data  in- 
eludes  failure  and  suspension  times  and  the  associated  theoretical 
distributions  to  be  assumed. 

The  fortran  listings  of  the  commuter  programs  are  lengthv; 

DISTSEL  and  RELIAB  contain  17  and  25  subroutines  respectively.  These 
listings  are  consequently  not  contained  in  this  report  but  can  be 
made  available  upon  request.  The  input  data  format  for  both  programs 
is  presented  in  Appendix  I  and  output  data  for  a  sample  problem  are 
given  in  Appendix  II. 

As  a  final  note,  it  should  be  emphasized  that  the  computer  models 
are  flexible.  Input  and  output  data  formats  can  readily  be  changed 
for  particular  situations,  perhaps  to  render  the  programs  compatible 
with  a  given  computer  data  file  svstem.  Improved  versions  of  the 
programs  will  also  be  generated  in  the  future. 
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APPENDIX  I 


INPUT  FORMAT  FOR  PROGRAMS  DISTSEL  AND  RELIAB 


The  following  listings  give  the  card  by  card  description  of  the 
innut  data  and  format  required  by  the  reliability  data  analysis 
programs . 


(a)  Overall  System  And  Subsystem  Data  Cards 


Card  Column 

No.  No.  Description  Format 


1 


1-80  Overall  system  identification. 


Alohanumeri c 


2  1-10 
11-20 
21-30 


3  1-40 


4  to  4  1-80 

+NSIJB 


Number  of  missions  over  system  life.  1 10 

Specified  system  life.  F10.0 

Desired  confidence  level  for  relia-  F10.0 

bility. 

Number  of  components  per  subsystem  4110 

for  which  data  is  supplied.  Up  to 
four  subsystems  can  be  handled. 


Subsystem  identification.  NSIIB  is  the 
number  of  subsystems  where  one  card 

is  used  to  describe  each  subsystem.  Alphanumeric 


(b)  Element  or  Component  Data  Cards 


The  following  data  cards  are  required  for  each  element 
starting  with  the  elements  of  the  first  subsystem  down  to  the  elements 
of  the  last  subsystem: 
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Card 

No. 


Column 

No. 


Description 


Format 


1 


1 

2-10 

11-80 


Theoretical  distribution  code  number*. 
Element  identification  code  number. 
Description  of  element. 


II 

19 

Alohanumeri c 


2  to  n^ 

♦  n„ 


2-11 

12-80 


Code  number  for  determining  if  data 
on  card  is  a  failure  time,  a  suspen¬ 
sion  time  or  last  card  for  current 
element**,  n ^  is  the  number  of 

failures  and  n  is  the  number  of 
s 

suspensions. 

Failure  or  suspension  time.  If  this 
is  the  last  card  then  the  remaining 
columns  are  blank. 

Description  or  the  failure  or  suspen¬ 
sion. 


II 


nl 0  %  0 


Alnhanumer^  r 


*  Distribution  Code  Number: 

1  -  Exponential 

2  -  Two  parameter  Weibull 

3  -  Three  parameter  Weibull 

4  -  Two  parameter  lognormal 

5  -  Three  parameter  lognormal 

6  -  Two  parameter  gamma 

7  -  Three  parameter  gamma 

8  -  Normal 

**  Failure  or  Suspension  Code  Number: 

1  -  Suspension 

2  -  Failure 

3  -  End  of  data 


(c)  Sample  Problem 

The  following  is  a  listing  of  data  input  for  a  hypothetical  reli¬ 
ability  data  analysis  problem: 
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SAMPLE  PROBLEM  TO  DEMONSTRATE 
10  2000.0  0.70 
3  1 

PROBLEM  SUBSYSTEM  A 
PROBLEM  SUBSYSTEM  B 
PROBLEM  SUBSYSTEM  C 


RELIABILITY  DATA  ANALYSIS  *'ODEL 


1 


SAMPLE 
SAMPLE 
SAMPLE 
1100000001 
2  1200.0 
000.0 
1500.0 
1300.0 
1350.0 
1000.0 
1300.0 
1100.0 
850.0 


2 
2 
3 

1100000002 

2  1500.0 
1  8000.0 

3 

1100000003 

l  onoo.o 

3 

2200000001 


SAMPLE  PROBLEM  ELEMENT  A1 
SPACE  FOR  FAILURE  DESCRIPTION 
SPACE  FOR  FAILURE 
SPACE  FOR  FAILURE 
SPACE  FOR  FAILURE 
SPACE  FOR  FAILURE 
SPACE  fqt?  faI LURE 
SPACE  FOR  FAILURE 
SPACE  FOR  FAILURE 
SPACE  EOF  FMLIPE 


DESCRIPTION 

DESCRIPTION 

DESCRIPTION 

DESCRIPTION 

DESCRIPTION 

DESCRIPTION 

DESCRIPTION 

DESCRIPTION 


SAMPLE  PROBLEM  ELEMENT  A7 
SPACE  fqr  FAILURE  DESCRIPTION 
SPACE  FOR  SUSPENSION  DESCRIPTION 

SAMPLE  PROBLEM  ELEMENT  A3 
SPACE  FOR  SUSPENSION  DESCRIPTION 

SAMPLE  PROBLEM  ELE'fENT  B1 


2 

1 200 . 0 

SPACE 

r0R 

rAILIJRE 

DESCRIPTION 

2 

ono.o 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

2 

1500.0 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

2 

1300.0 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

2 

1350.0 

SPACE 

^OR 

FAILURE 

DESCRIPTION 

2 

1000.0 

SPACE 

nQR 

FAILURE 

DESCRIPTION 

2 

1300.0 

SPACE 

POR 

FAILURE 

DESCRIPTION 

2 

1100.0 

SPACE 

FOR 

failure 

DESCRIPTION 

2 

850.0 

SPACE 

FOP 

FAILURE 

DESCRIPTION 

4300000001 

SAMPLE  PROBLEM  ELEMENT  Cl 

2 

1200.0 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

2 

000.0 

SPACE 

FOR 

failure 

DESCRIPTION 

2 

1500.0 

SPACE 

r0R 

failure 

DESCRIPTION 

2 

1300.0 

SPACE 

FOR 

failure 

DESCRIPTION 

2 

1350.0 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

2 

1000.0 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

2 

1300.0 

SPACE 

tor 

FAILURE 

DESCRIPTION 

2 

1100.0 

SPACE 

nOR 

FAILURE 

DESCRIPTION 

2 

850.0 

SPACE 

non 

FAILURE 

DESCRIPTION 

5P 


APPENDIX  II 


OUTPUT  DATA  TT)R  SAMPLE  PROBLEM 


The  following  pages  contain  outnut  data  ronerated  hv  nronrams 
DISTSEL  and  RELIAB  for  the  sannle  input  data  presented  in  Appendix  I. 
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RESULTS  FROM  PROORAM  DISTSEL 
FAILED  SPECIMENS  AND  MEDIAN  RANKS.,. 


850.0 

0,07447 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

900,0 

0.18085 

SPACE 

HOP 

FAILURE 

DESCRIPTION 

1000.0 

0.28723 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

1100.0 

0.39362 

SPACE 

PO" 

FAILURE 

DESCRIPTION 

1200.0 

0,50000 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

1300.0 

0.60638 

SPACE 

POR 

FAILURE 

DESCRIPTION 

1300.0 

0.71277 

SPACE 

POR 

FAILURE 

DESCRIPTION 

1350.0 

0.81915 

SPACE 

FOR 

FAILURE 

DESCRIPTION 

1500.0 

0.92553 

SPACE 

POR 

FAILURE 

DESCRIPTION 
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STANDARD  ERROR  0^  ESTP1ATE . . . 


DISTRIBUTION  S.E. 


EXPONENTIAL 
NORMAL 
LOO-NORMAL 
WEI  BULL-2  PARM 
WEI  BULL- 3  PARM 
GAMMA- 2  PARM. 


0.280808E  00 
0.566338E-01 
0.657507E-01 
0.530806E-01 
0.64220SE-D1 
D.832000E-01 


DISTRIBUTION  SELECTED...  WEIBULL-2  PARAMETER 

K-S  STATISTIC  (  0.141)  IS  LESS  THAN  TABLE  VALUE  (  0.388) 

THEREFORE  MAY  ACCEPT  HYPOTHETICAL  DISTRIBUTION  AT  10  PERCENT  LEVEL 
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SAMPLE  PROBLEM  TO  DEMONSTRATE  RELIABILITY  ORTA  ANALYSIS  MODEL 


MARCH  1973 


0DRTR 


lOOOOOOOO  SAMPLE  PROBLEM  ELEMENT  fll 

EXPONENTIAL  D ISTR IBUT ION. .. LEAST  SOUAHES  ESTIMATE  MEAN  LIFE  =  HUB. 173 


*3 


CUMULATIVE  DISTRIBUTION 


SAMPLE  PROBLEM  TO  DEMONSTRATE  RELIABILITY  DATA  ANALYSI 5  MODEL 


MARCH  1973 


DATA 

1500 


1370 


1240 


1110 


980 


850 


100000000  SAMPLE  PROBLEM  ELEMENT  Al 

NORMAL  DISTRIBUTION. . .LEAST  SQUARES  ESTIMATES  MEAN  LIFE 
STANDARD  DEVIATION  =245.0730 


1166.665 


SAMPLE  PROBLEM  TO  DEMONSTRATE  RELIABILITY  DATA  ANALYSIS  MODEL 


MARCH  1973 


DATA 

100000000  SAMPLE  PROBLEM  ELEMENT  R1 

LOG  NORMAL  01 STRI BUT  1  ON. . . LEAST  SQUARES  ESTIMATES  MEAN  LIFE  -  1175.337 

STANDARD  DEVIATION  =259.1138 

1500  _ 

r  /  x 

1339  _ 

/  K 

X  X 

1195  _ 

r  X  / 

1067  _ 

952 

- 1 

X 

aso 

X  ' 

.001  .01  .05  .10  .20  .10  .60  .BO  .90  .95  .99 

CUMULATIVE  DISTRIBUTION 

.  US  .10 
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DRTfl 

SAMPLE  PROBLEM  TO  DEMONSTRATE  RELIABILITY  DATA  ANALYSIS  MODEL  MARCH  1973 

1U0D0DD0D  SAMPLE  PROBLEM  ELEMENT  A1 

2-PARM.  HEIBULI _ LEAST  SQUARES  ESTIMATES  SHAPE  =  5.587 

SCAi_E  =  1259.  OD 

1500  _| 

f 

j 

1339  _ 

/  X 

X  /x 

1195  - 

J 

1067  _ 

952  _ 

j  K 

/  x 

850  _ 

X 

.001  .01  .05  .10  .20  .40  .60  .80  .90 

CUMULATIVE  DISTRIBUTION 


.01 
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sample  problem  to  demonstrate  reliability  DATA  ANALYSIS  model 


MARCH  1973 


DATA 

903 


700 


543 


421 


326 


253 


100000000  SAMPLE  PROBLEM  ELEMENT  A1 

3-PARM.  NEIBULL... LEAST  SOUAHES  ESTIMATES  SHAPE  =  2.494 
SCALE  =  648. B66 


.001  .01  .05  .10  .20  .40  .60  .00  .90 

CUMULATIVE  DISTRIBUTION 
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CUMULRT I  V  E  DISTRIBUTION 

.00  0.13  0.25  0.36  0.50  0.63  0.75  0.88  1.00 
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RESULTS  OF  COMPUTATIONS  FOR  ELE'tENT  ID  NIP1BER  100000001 

SAMPLE  PROBLEM  ELEMENT  A1 
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STANDARD  DEVIATION  =  1166.6666666667 
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RESULTS  OF  COMPUTATIONS  POR  ELEMENT  ID  N1P1BER  100000002 

SAMPLE  PROBLEM  ELEMENT  A2 
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RESULTS  OF  COMPUTATIONS  FOR  E LENIENT  ID  NIHBER  100000003 

SAMPLE  PROBLEM  ELEMENT  A3 
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SCALE  PARAMETER  =  1253.3223250387 

SHAPE  PARAMETER  =  6.5319380785 

AVERAGE  =  1168.1272794997 

STANDARD  DEVIATION  =  209.2752304140 
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INTERVAL  RELIABILITY  FOR  EVERY  MISSION 

0.9999927  0.9994249  0.9924625  0.9562320  0.8474288 

0.6753514  0.6540568  0.8713890  0.9602432  0.91183*6  (Continued) 


COMPONENT  CONFIDENCE  LIMIT  ON  RELIABILITY  IS  PRESENTLY 
DETERMINED  ASSUMING  CONSTANT  FAILURE  RATE 
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RESULTS  OP  COMPUTATIONS  POR  ELEMENT  ID  NUMBER  .300000001 
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SCALE  PARAMETER  *  7.D45S385748 

SHAPE  PARAMETER  =  0.1829129466 

AVERAGE  =  1167.0881299857 

STANDARD  DEVIATION  =  215.2736078302 
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COMPONENT  CONFIDENCE  LIMIT  ON  RELIABILITY  IS  PRESENTLY 
DETERMINED  ASSUMING  CONSTANT  FAILURE  RATE 
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SUMMARY  OF  ELEMENT  RELIABILITY  COMPUTATIONS 
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SUBSYSTEM. . .SAMPLE  PROBLEM  SUBSYSTEM  A 

mooonnoi  sample  problem  element  ai 

100000002  SAMPLE  PROBLEM  ELEMENT  A2 

10000000.3  SAMPLE  PROBLEM  ELE'fENT  A3 


SUBSYSTEM...  SA‘1PLE  PROBLEM  SUBSYSTEM  B 

200000001  SA»1PLE  PROBLEM  ELEMENT  B1 


I 


CO 


a 


83 


VC 

to 

Tf 

rr 

to 

to 

d 

c 

d 

vC 

c 

d 

c 

rr 

oc 

vC 

»— < 

c 

C- 

LO 

rr 

c 

u 

e 

to 

cc 

c 

c 

oc 

vC 

c 

c 

rr 

cr 

oc 

LT. 

d 

r 

d 

oc 

c 

oc 

oc 

>- 

d 

» 

>- 

oc 

oc 

>- 

CO 

* 

> 

• 

c 

• 

» 

u 

T. 

c 

t 

c 

c 

tr 

cr  c 

c 

H 

C 

c 

c 

HH 

UJ 

cr 

UJ 

cr 

U- 

H 

d 

UJ 

c 

It 

d 

UJ 

e 

or 

UJ 

<f 

II 

d 

e 

II 

d 

C  It 

II 

d 

E- 

t 

|7 

<f 

o 

l  • 

< 

cr 

< 

C_J 

Cl 

c 

> 

C_) 

d 

c 

CJ 

d 

a. 

>- 

c: 

>* 

H 

££ 

>- 

T-' 

h- 

CO 

UJ 

cr 

H 

KH 

to 

UJ 

p 

uc 

UJ 

C 

HH 

HH 

Cl 

HH 

HH 

d 

HH 

HH 

d 

C-J 

UJ 

• 

UJ 

►H 

•  u 

UJ 

r 

hH 

c 

HH 

cc 

>- 

b 

C  HH 

HH 

>- 

>* 

C3 

E- 

HH 

CO 

< 

H 

hH 

CO 

CO 

E— 

H-< 

h- 

< 

HH 

< 

< 

►H 

►H 

< 

< 

< 

HH 

< 

HH 

i-h  HH 

UJ 

u 

HH 

CM 

UJ 

UJ 

u 

HH 

to 

HH 

UJ 

u 

UJ 

-J 

d 

ii 

UJ 

UJ 

l-H 

II  UJ 

UJ 

HH 

HH 

a1  uj 

cc 

H 

UJ 

d 

d 

CO 

H 

UJ 

d 

UJ 

CC 

OQ 

uj  d 

< 

d 

UJ 

< 

d 

UJ 

Cr 

< 

r' 

< 

c: 

►H 

-J 

m 

►H 

UJ 

cc 

HH 

< 

►H 

UJ 

w 

f* 

UJ 

H 

uj  r: 

T. 

UJ 

E— * 

-J 

b  o 

UJ 

co 

> 

o 

c 

UJ 

co 

>  o 

C 

UJ 

CO 

UJ 

2:  »-h 

d 

UJ 

►H 

i 

d 

r 

UJ  HH 

HH 

d 

rr 

Cl 

CO 

o 

*-J 

CO 

CO 

o 

^-J  CO 

to 

c 

r  co 

UJ  t-H 

Q 

CJ 

CO 

to 

cr. 

C_) 

to 

T” 

to 

cr 

u 

UJ 

UJ 

HH 

UJ 

HH 

u 

UJ  tH 

UJ 

HH 

UJ 

UJ 

h  r 

C_) 

cr 

o 

0f 

h 

Jr 

u 

cr 

u  r 

H 

r ' 

u 

c 

fc 

co 

CO 

n 

2: 

PO 

to 

rr 

>-  u 

UJ 

HH 

UJ 

a 

>- 

U 

UJ 

►H 

UJ  c 

>- 

u 

UJ 

HH 

>- 

CO  c 

c. 

cr 

uj 

CO 

c 

cr 

J~ 

r  uj 

to 

c 

r 

CO 

03 

HH 

P 

HH 

u 

oc 

HH 

b 

u 

CO 

HH 

OQ 

rr  w 

U 

CO 

u 

sr 

rr 

UJ 

U 

co 

u  n 

rr 

UJ 

U 

to 

p 

CO  E- 

2: 

CO 

r: 

UJ 

CO 

E- 

CO 

n  uj 

to 

H 

rr 

to 

CO 

< 

o 

< 

o 

C 

< 

C 

< 

c  r 

< 

c 

< 

Cl  T 

O 

u 

HH 

d 

CJ 

cu  ^ 

d 

T~ 

cu 

UL 

O  hh 

cr 

u 

o 

KH 

r 

U- 

o 

HH 

cr 

o 

U.  E- 

y 

UJ 

e 

2: 

u 

H 

UJ 

£T  n 

u 

E— ' 

UJ 

CO 

UJ 

UJ 

c 

to 

UJ 

UJ  o 

CO 

UJ 

CO 

>-  UJ 

H 

hh 

HH 

u 

>~ 

UJ 

H 

hH 

HH  c_j 

>- 

UJ 

E— 

HH 

H 

E— » 

CO 

y 

U. 

to 

J“ 

U 

E- 

to 

Hr 

UJ 

►H  d 

>- 

Cl 

HH 

d 

WH 

H 

>- 

d 

hh  d 

HH 

>- 

d 

J  ' 

uj  s 

CO 

UJ 

u 

to 

UJ 

CJ  UJ 

-J 

to 

UJ 

CO 

HH  HH 

OQ 

H 

UJ 

hH 

►— i 

OQ 

fe 

UJ  3r 

HH 

HH 

CO 

E- 

UJ 

CQ  O 

UJ 

Dm 

o 

CO 

o 

rr 

d  o 

GQ 

o 

r 

UJ 

d 

<  d 

CO 

a 

CO 

hJ 

< 

a. 

CO 

c 

CO  UJ 

< 

d 

to 

c 

uj 

d 

o 

UJ 

u 

UJ 

Q 

►H 

U. 

o 

u 


UJ 

UJ 

d 

C 

UJ 

U 

UJ 

c 

HH 

u, 

o 

u 


a 

C: 

UJ 

LJ 

y. 

uj 

c 

HH 

u. 

o 

LJ 


84 


SPECIFIED  CONFIDENCE  LEVEL  =  0.70000000^" 

LOWER  CONFIDENCE!)  MISSION  RELIABILITY  =  n.90S0]ni007 
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